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The antineutrino spectra measured in recent experiments at reactors are inconsistent with calcu¬ 
lations based on the conversion of integral beta spectra recorded at the ILL reactor. ®^Rb makes 
the dominant contribution to the reactor antineutrino spectrum in the 5-8 MeV range but its decay 
properties are in question. We have studied ®^Rb decay with total absorption spectroscopy. Previ¬ 
ously unobserved beta feeding was seen in the 4.5-5.5 region and the GS to GS feeding was found 
to be 87.5(25)%. The impact on the reactor antineutrino spectra calculated with the summation 
method is shown and discussed. 


Beta decay properties of fission products are at the 
origin of the antineutrino flux emitted by reactor cores. 
This flux has been used for decades as a source for reactor 
neutrino experiments, such as Daya Bay, Double Ghooz 
and Reno which have recently published their new results 
for the mixing angle 0 i 3 [IH3]. These results will allow 
future searches for the CP violation phase 5 or the neu¬ 
trino mass hierarchy with complementary experiments 
at reactors @] . The accurate determination and under¬ 
standing of the emitted reactor antineutrino flux is thus 
still required for present and future experiments. The 
recent re-estimate of reactor antineutrino energy spec¬ 
tra [SJ [5] has led to the so-called “reactor anomaly” [7] , 
at the origin of new experimental projects chasing short 
distance oscillations at research reactors [8] . These calcu¬ 
lations are based on the conversion into antineutrinos of 
the only available measurement of the beta energy spec¬ 
tra performed using the high flux reactor at the Insti- 
tut Laue-Langevin (ILL) in Grenoble, France |1]. The 
conversion method has until now been considered as the 
most precise one by experimenters studying neutrino os¬ 


cillations. But recently, Hayes et al. m have shown that 
it is dependent on the underlying nuclear physics and 
that the associated errors should be revised. In addition, 
very recent experimental results from m have shown 
an unexplained distortion between 4 and 8 MeV in their 
measured positron energy spectra from Pressurized Wa¬ 
ter Reactors (PWR) |TT] with respect to the converted 
spectra EE] (an excess of ca. 10% over 2 MeV followed 
by a dip). The positron energy used in reactor antineu¬ 
trino experiments corresponds to the antineutrino energy 
minus the mass difference between the neutron and pro¬ 
ton. In this context, new evaluations of PWR antineu¬ 
trino energy spectrum are essential. 

An alternative method, independent of the ILL mea¬ 
surements, relies on the summation of the contributions 
of the fission product beta decay branches to obtain the 
antineutrino energy spectra. The need to measure new 
nuclear physics properties of some major contributors to 
the antineutrino spectrum was underlined in |12j , where 
it was shown that we should use the Total Absorption 
Spectroscopy (TAS) technique to avoid the pandemo- 
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nium effect m and improve the predictions of the sum¬ 
mation method. The summation method is indeed the 
only one which allows the prediction of antineutrino spec¬ 
tra for which no integral beta measurement exists. This 
is required, for instance, in the context of the R&D of an¬ 
tineutrino detection as a tool for reactor monitoring El- 

In this Letter, we present the hrst results of an exper¬ 
imental campaign performed with TAS technique [15] 
aimed at the measurement of beta decay properties of 
important contributors to the reactor antineutrino en¬ 
ergy spectrum emitted by PWRs. In particular, we show 
new results for ®^Rb, the largest contributor to the reac¬ 
tor antineutrino flux in the energy range above 5 MeV. In 
the following, we present a short list of nuclei making the 
main contributions to the antineutrino energy spectrum 
above 4 MeV, obtained using the summation method pre¬ 
sented in m- Then, previous experimental knowledge of 
the beta decay properties of ®^Rb is summarized, and the 
TAS method, the experimental setup used and the data 
analysis performed are presented. Finally, we show the 
beta feeding obtained and present the impact of the new 
results on reactor antineutrino energy spectra. 

The main contributors to the antineutrino energy spec¬ 
tra from 4 to 8 MeV are listed in TablelU In our calcula¬ 
tion, we have chosen to minimize the impact of the pan¬ 
demonium effect on the antineutrino spectra and on the 
computed proportions of the nuclei per energy bin. For 
this purpose, we have used data from [TS] for 

^^^Cs, ^^®Te and from [T7j for ®®Sr and ®°Rb, be¬ 
cause these two sets of data are likely to be pandemo¬ 
nium free (though they may suffer from other systematic 
errors). As they were not measured by m or m. data 
for ssmy lo^Nb were taken from [TS] and 
from m- Indeed, a careful choice of data sets is needed, 
especially to select nuclei which would deserve new mea¬ 
surements, as is illustrated below with the case of ®^Rb. 
Note that in the 4 to 6 MeV range, unknown nuclei re¬ 
quiring the use of models represent less than 1% of the 
spectrum, while they represent about 4% of the 7 to 
8 MeV bin. 

®^Rb makes the main contribution between 4 to 8 MeV, 
representing alone up to about 38% of the 7 to 8 MeV 
bin and 16 % of the 5 to 8 MeV range. ®^Rb is quite 
controversial: the beta feeding to the ground state of its 
daughter nucleus, ®^Sr, was fixed at 51% ± 18% in the 
ENSDF data base [20] until 2012, before the inclusion in 
the references of the article from Lhersonneau et al. |2T| , 
which concluded that close to half of the decay intensity, 
mostly high energy ground state transitions, is missing 
in the decay scheme. Following this reference the beta 
feeding to the ground state of ®^Sr was recently changed 
to 95.2 % ± 0.7% in the ENSDF database [l9|. 

®^Rb has a large <5/3 value which makes it a good 
candidate to be a pandemonium nucleus. The pande¬ 
monium effect m arises from the difhculty encountered 
in building level schemes for complex beta decays using 


Germanium detectors, especially when beta transitions 
occur to high-energy levels or regions of high level den¬ 
sity. This leads to an underestimate of the correspond¬ 
ing beta branches to states at high excitation energy and 
thus to a distortion in the beta decay feeding. In addi¬ 
tion, ®^Rb has also been used as a critical example [55] 
to show how beta-decay strength calculations impact on 
the predictive power of models in reconstructing half- 
lives and beta-delayed neutron emission probabilities of 
nuclei, whose properties are important in the simulation 
of the astrophysical r-process. It is also on NEA/IAEA 
lists of important contributors to reactor decay heat [33] . 

A total absorption spectrometer is a calorimeter mea¬ 
suring the gamma cascades emitted by the deexcitation 
of the daughter nucleus after beta decay of the parent. 
The detection of the total energy allows the deduction of 
the feeding probability of excited levels populated in the 
beta decay. This quantity is calculated by solving the 
“inverse problem” as presented below. The beta feeding 
/ gives direct access to the beta intensity li = fi/T^kfk 
and then to the beta strength, a microscopic quantity 
that can be directly compared with models [24]. The 
detector used in the measurement of the ®^Rb decay is 
composed of 12 crystals of BaF 2 arranged in a compact 
geometry described in [3^. Each crystal is coupled to 
a photomultiplier tube converting the scintillation light 
into an electrical signal directly proportional to the de¬ 
tected energy. The gamma detection efficiency was ^ 
80 % at 5 MeV. This spectrometer was coupled to a sili¬ 
con detector placed in the center, behind the source im¬ 
plantation zone, to tag the beta emission. This reduces 
the background by demanding coincidences between beta 
events and the following gamma emission from deexcita¬ 
tion of levels in the daughter. 


TABLE I. Main contributors to a standard PWR antineu¬ 
trino energy spectrum computed with the MURE code cou¬ 
pled with the list of nuclear data given in assuming that 
they have been emitted by ^^^U (52%), ^®®Pu (33%), ^^^Pu 
( 6 %)and (3 7 %) for a 450 day irradiation time and using 

the summation method described in | 12 | . 



4 - 5 MeV 

5 - 6 MeV 

6 - 7MeV 

7 - 8 MeV 

®^Rb 

4.74% 

11.49% 

24.27% 

37.98% 

96y 

5.56% 

10.75% 

14.10% 

- 

i^^Cs 

3.35% 

6 . 02 % 

7.93% 

3.52% 

loONb 

5.52% 

6.03% 

- 

- 

®3Rb 

2.34% 

4.17% 

6.78% 

4.21% 

98rn, Y 

2.43% 

3.16% 

4.57% 

4.95% 

i35Te 

4.01% 

3.58% 

- 

- 

104mNb 

0.72% 

1.82% 

4.15% 

7.76% 

®°Rb 

1.90% 

2.59% 

1.40% 

- 

^®Sr 

2.65% 

2.96% 

- 

- 

®"Rb 

1.32% 

2.06% 

2.84% 

3.96% 
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ions were produced via proton-induced fission on 
an uranium target at the IGISOL facility [26] in the ac¬ 
celerator Laboratory of the University of Jyvaskyla (Fin¬ 
land). JYFLTRAP double Penning trap |17| was used 
for selecting with high precision only ®^Rb ions using the 
mass-selective buffer gas cooling technique |2H]- This 
high level of purification of the beam is necessary in TAS 
experiments in order to reduce systematic uncertainties 
related to the purity of the beam. 

As stated above the main observable in a TAS mea¬ 
surement is the beta feeding to the energy levels of the 
daughter nucleus which is contained in the measured 
gamma spectrum convoluted with the detector response. 
To extract this information we have to solve the so-called 
“inverse problem”. It consists of solving the equation 
di = ' X /j, where is the response matrix of the 

detector to an assumed decay level scheme {B). Rij con¬ 
nects feeding to level j (/j) to counts in the bin i of the 
“measured” TAS spectrum (di). The analysis procedure 
has been described in previous publications [Ml HH] and 
is well understood. 

To perform the analysis, very clean decay data d are 
required. Possible contamination from the daughter nu¬ 
clei decay and pileup signals are subtracted from the raw 
data. The shape of the pile-up spectrum has been com¬ 
puted by summing in the ADC time window two events 
randomly extracted from the raw data. The absolute 
normalization of the pileup was performed using the data 
counting rate and the ADC time window in which ran¬ 
dom coincidences can occur [301 in]. The shape of the 
spectrum from the decay of the daughter nucleus ®^Sr 
has been simulated using its known level scheme from 
ENSDF [Tni and the detector response. The normal¬ 
ization factor has been obtained solving the Bateman 
equations for the decay of ®^Rb in realistic experimen¬ 
tal conditions, i.e. considering the experimental time for 
implantation and measuring cycles. The contamination 
from ®^Sr decay represents 0.08 % of the total ®^Rb data 
acquired. The response matrix R is calculated by simu¬ 
lating the detector response to beta and gamma cascades 
emitted during the decay with a dedicated CEANT4 [31] 
Monte Carlo simulation. The latter has been validated 
using measurements performed with known sources in 
order to reproduce the detector response in great detail 
[3011351 133] ■ The “inverse problem” is solved by using a 
maximization expectation algorithm based on the Bayes 
theorem and combined with a minimization [35] . It 
makes use of an iterative method to find the final feeding 
distribution by minimizing the difference between the ex¬ 
perimental data and the spectrum recreated by the result 
of the algorithm at each iteration. The analysis starts 
with a first guess at feeding values extracted from the 
literature, or an equally probable feeding distribution if 
the nucleus is poorly known, and stops when the value 
deduced from the two spectra is at a minimum. 

The starting point for solving the inverse problem is 


the construction of the branching ratio matrix (B) for the 
states populated in the decay. For this purpose we begin 
by using the known information, derived from high reso¬ 
lution studies, about levels up to an excitation energy of 
1778 keV in ®^Sr [T3]. Above this energy little is known 
and the data are divided into 40 keV bins up to the Qp 
value. In this range we must have recourse to semiem- 
pirical statistical models and we must supply as input 
both the level densities and gamma strength functions. 
Three level-density models were tested: Back-Shifted- 
Fermi-Cas (BSFC) [3HI3Z], Constant-Temperature [37] . 
and Cilbert-Cameron models [31]. The last of these is 
a combination of the other two. The Cilbert-Cameron 
formulation was chosen because it best reproduces the 
experimental data at low energies. The gamma strengths 
were modeled with a Lorentz function using the param¬ 
eters given in [33]. In determining the /3-feeding distri¬ 
bution, it is possible to fix or vary the feeding to each 
individual level or energy bin. The feeding to the 1673.3 
keV level was set to zero, since the probable spin parity 
is 4+ and any feeding from the ®^Rb ground state must 
be negligible. 

The reconstructed spectrum (blue dashed line) calcu¬ 
lated using the feeding distribution obtained from this 
analysis is compared with the clean decay data (black 
continuous line) of ®^Rb in the upper panel of figure[^ 
[30] . The lower panel shows the residues between these 
two curves. The beta intensity obtained from the so¬ 
lution of the inverse problem for ®^Rb is shown in fig¬ 
ure]^ in the blue continuous line, while the red dashed 
lines are the intensities from ENSDF [19]. As the ground 
state feeding is very important in the case of the decay 
of ®^Rb, we have estimated the main errors involved in 
this reconstruction [30]. They are listed as follows: the 
threshold of the beta spectrum, statistical uncertainty, 
error induced by pile-up subtraction, errors in the detec¬ 
tor energy calibration and resolution used in the calcu¬ 
lation of the response matrix R and errors obtained by 
testing different input parameters for the calculation of 
R and inverse problem resolution. A sum in quadrature 
of all the systematic and statistical errors quoted above 
gives a 2.5 % error on the ground state feeding. This re¬ 
sult is conservative, as we have voluntarily adopted large 
values of the main errors which are associated with the 
threshold of the beta spectrum and with the choice of 
model for the level density. The TAS results show some 
beta intensity around 4.5 and 5.5 MeV which was not 
detected before. The intensity to the ground state ob¬ 
tained from our analysis is 87.5 (25) %. This value can 
be obtained from the data analysis because the TAS de¬ 
tector also measures the bremsstrahlung radiation from 
the beta particles. These events are in the low energy 
part of the measured spectrum and, since they are con¬ 
sidered in the response matrix R, they contribute to the 
reconstruction of the spectrum and, then, in the calcula¬ 
tion of beta feeding. The selected ground state feeding is 
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FIG. 1. Upper panel: Comparison between measured spec¬ 
trum (black continuous line) and reconstructed one (blue 
dashed line) with the feeding obtained from the TAS data 
analysis. Lower panel: Residues between the two curves re¬ 
ported in the upper panel. 



FIG. 2. Beta intensity for the decay of ®^Rb obtained with 
TAS measurements in the blue continuous line. Red dashed 
lines are data from ENSDF [19]. Ground state feeding from 
ENSDF is 95.2% ± 0.7% , while from this work it is 87.5% ± 
2.5% 


the one which minimizes the value determined from 
the experimental data and the reconstructed spectrum 
after the analysis. If we fix the ground state feeding to 
be 95.2 % as reported in the ENSDF data base [TD our 
analysis converges with a value of 2048 which is much 
larger than the minimum of 630, completely excluding 
this hypothesis. 

The antineutrino energy spectrum emitted in ®^Rb 
beta decay has been computed using the beta feeding pre¬ 
sented above. The GS to GS transition is first-forbidden 
nonunique (spin parity of ®^Rb 0“). Different spectral 
shapes were assumed for this transition, considering the 
different possibilities listed in uni; an allowed shape, a 
first forbidden nonunique shape due to the GT operator 


and a first forbidden nonunique shape due to the pA oper¬ 
ator. It was also assumed that the remaining transitions 
were of allowed or first forbidden unique type. The vari¬ 
ous combinations of these options were computed and the 
shapes obtained were very similar. No significant impact 
is expected from the uncertainty of the shape of the first 
forbidden nonunique GS to GS transition. We chose to 
adopt an allowed shape for the GS to GS transition and 
first forbidden unique shapes for the remaining branches 
due to the spins and parities of the known transitions in 
this nucleus. 

The antineutrino energy spectra were calculated with 
the summation method described in |12j . In m the data 
adopted for ®^Rb were extracted from [16] . In principle, 
these measurements should not suffer from the pandemo¬ 
nium effect, nor from a lack of knowledge of the types of 
the beta transitions. Unfortunately, however, the error 
bars are quite large. In figurej^ the ratio between the an¬ 
tineutrino spectra of 239,24ip^ 235,238pj fj-Qjjg [^2] and 

those obtained using our new results for ®^Rb is displayed 
with the red dashed-dotted line. As expected, the main 
effect is in the 4 to 8MeV antineutrino energy range, 
with a maximum between 7 and 8 MeV, and amounts to 
4.5% for 235 u, 3.5% for ^sspu, 2% for 24ipu and 1.5% for 
238u These discrepancies are due to the difference in the 
shapes of the antineutrino spectra built with the newly 
measured beta feedings with respect to the antineutrino 
spectra converted from Rudstam’s measurements. The 
comparison would be very similar if we had used the lat¬ 
est ENSDF [19] data for ®^Rb in our summation calcula¬ 
tions, as was done in [40]. The ratio is displayed as well in 
Fig.|3] with green dotted lines, and is nearly superposed 
on the ratio built when using Rudstam data in the first 
place. The change becomes even more dramatic if one 
compares with summation method spectra in which an 
older version of the ENSDF data was used, as in [41]. The 
latter ratio is plotted with black dashed lines in figurej^ 
This shows the relevance of the present ®^Rb decay data 
in the calculations. 

In summary, the results of new measurements of the 
beta decay properties of ®^Rb have been presented. 
This nucleus makes one of the largest contributions 
to the emitted antineutrino flux by standard thermal 
reactors in the energy region above 5 MeV. The mea¬ 
surements have been performed using pure isotopic 
beams and the TAS technique to provide data free 
from the pandemonium effect. The measured feeding 
distribution, which extends to states previously not seen 
in high-resolution measurements and also determines the 
GS to GS feeeding as 87.5(25)%, confirms the relevance 
of this decay to antineutrino summation calculations. 
The impact of the measurements has been evaluated 
by comparing the ratio of summation calculations using 
the new feeding distribution with the results using the 
feeding distributions employed in [H] SOI m] . The effect 
of introducing the new results is particularly marked in 
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Energy (MeV) 

FIG. 3. Ratio between the antineutrino spectra calculated 
using the results presented in this paper with respect to the 
data on decay used in [T^ (thick red dashed-dotted line), 
in |40| (green dotted line) and in |41| (black dashed line). The 
sharp drop in the ratio, in one single bin located at the Q 
value of the ®^Rb, is due to the different values in Q given 
in m and [19] , that were used to reconstruct the antineu¬ 
trino spectrum. A gray horizontal bar is placed above the 
antineutrino energy scale to indicate the region of the distor¬ 
tion observed by the reactor antineutrino experiments whit 
respect to converted spectra. 


the case of [41] and calls for a revision of the conclusions 
drawn in that paper. It is clear that this is because 
the GS to GS feeding used in [IT] was incorrect. The 
overall agreement of the new summation calculations 
with the converted spectra |6] is improved in the 4 
to 8 MeV range except in the case of ^35^ j-pj. 
the summation method spectrum is always below the 
converted spectrum. The change is especially striking in 
the case of m in the 5 to 8 MeV antineutrino energy 
range, which overlaps the energy region in which reactor 
neutrino experiments have shown a spectral distor¬ 
tion m- It also shows that the inclusion of all existing 
TAS nuclear data (ca. 37 nuclei) in |?T] ’s calculation may 
change dramatically the spectral shape they compute. 
Overall, this emphasizes why new measurements are 
needed for the radioactive decays of importance in 
the reactor antineutrino spectrum and, in particular, 
why measurements should be performed with the total 
absorption method. The present measurement, which 
reduces significantly the uncertainties associated with 
the antineutrino summation calculations in the 4 to 
8 MeV range, is an important step towards better pre¬ 
dictions with the summation method. Provided that in 
the 4 to 6 MeV range, unknown nuclei requiring the use 
of models represent less than 1% of the spectrum in the 


summation calculation from nn, one can thus expect a 
dramatic reduction of the final uncertainty in this range 
as a long term result of the TAS campaign. In parallel, 
the impact of the uncertainties of the fission yields on 
the antineutrino spectrum needs to be evaluated more 
accurately. 
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